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The Sivers transverse single spin asymmetry (TSSA) is calculated and compared at different 
scales using the TMD evolution equations applied to previously existing extractions. We apply the 
Collins-Soper-Sterman (CSS) formalism, using the version recently developed by Collins. Our calcu- 
lations rely on the universality properties of TMD-functions that follow from the TMD-factorization 
theorem. Accordingly, the non-perturbative input is fixed by earlier experimental measurements, 
including both polarized semi-inclusive deep inelastic scattering (SIDIS) and unpolarized Drell-Yan 
(DY) scattering. It is shown that recent preliminary COMPASS measurements are consistent with 
the suppression prescribed by TMD evolution. 



Conventional collinear perturbative QCD (pQCD), 
when applied in its range of applicability, has proven suc- 
cessful for over three decades Along the way, it has 
illustrated the importance of evolution for relating physi- 
cal observables to fundamental quark-gluon QCD degrees 
of freedom in a unified formalism. The classic applica- 
tions of pQCD require parametrizations of collinear par- 
ton distribution functions (PDFs), wherein all intrinsic 
transverse motion of the confined partons is neglected 
inside the hard part of the collision and integrated over 
in the definitions of the PDFs. The PDFs contain in- 
formation about the intrinsic non-perturbative structure 
of the hadron and have clear operator definitions with 
well- understood scale dependence (QCD evolution). 

An important next step is to achieve a similarly 
successful application of perturbative QCD that takes 
into account the emerging picture of the hadron as a 
three-dimensional dynamical object composed of par- 
tons with their own intrinsic motion. In such studies, 
the relevant observables are not properly handled by 
standard collinear factorization, and one is confronted 
with the objects like TMD parton distribution func- 
tions (TMD PDFs) and TMD fragmentation functions 
(FFs). (Collectively, we refer to such objects as "TMDs" ) 
These, like their collinear counterparts, describe the non- 
perturbative properties of the external hadrons; but, un- 
like the collinear PDFs, the TMDs also account for the 
intrinsic transverse Fermi motion of the bound partons. 

Several experimental facilities, including HERMES 
(DESY), COMPASS (CERN) and JLab, explore these 
distributions. Moreover, the future 12 GcV JLab up- 
grade and a planned Electron Ion Collider (EIC) [2j 
will provide new opportunities to experimentally probe 
hadron structure. A particularly interesting TMD PDF 
is the so-called Sivers function, which is interpreted as 
the probability density for finding a parton with a given 
transverse and longitudinal momentum inside a trans- 
versely polarized hadron (usually proton) target. That 
it arises at leading power in pQCD is due to interesting 



non-perturbative aspects of QCD related to time rever- 
sal and parity invariance. In polarized SIDIS, it gives 
a sin^/j — (/>s) azimuthal modulation to the differential 
cross section, cf>s and 4>h being the azimuthal angles of the 
initial transverse hadron spin and the final state hadron 
transverse momentum respectively. This letter will focus 
on a comparison of recent theoretical treatments of the 
Sivers function with recently available experimental data 
on TSSAs in SIDIS, and predict the size of the asymme- 
try for future extractions at larger Q. 

Like the collinear PDFs, TMDs evolve with the hard 
scale Q. To properly account for this it is imperative to 
work in a QCD factorization formalism that incorporates 
well-defined TMDs. The original TMD-factorization for- 
malism was developed by Collins, Soper, and Sterman 0- 
5] in the context of e + e~-annihilation to back-to-back 
jets and the unpolarized Drell-Yan process. The CSS for- 
malism was later extended by Ji, Ma and Yuan to SIDIS 
in Ref. |(|, and to include polarization in Ref. 0]. The 
application of Collins-Soper (CS) evolution was extended 
to the spin-dependent case by Idilbi et al in Ref. [|| . Fur- 
thermore, an implementation of a CS-style evolution has 
been applied in Ref. [9] to the calculation of TSSAs by 
including a Sudakov form factor, leading to approximate 
power-law Q-behavior for the peak of the asymmetry. 

For quite some time, however, a satisfactory treatment 
of TMD-factorization remained incomplete because of a 
lack of good definitions for the TMD PDFs and FFs 
themselves The most commonly quoted defini- 

tions suffered from unphysical divergences, rendering it 
unclear which objects should be used in parametrizations 
of experimental data and treated as universal in the usual 
sense of a factorization theorem. This is particularly 
problematic for studies that purport to extract proper- 
ties intrinsic to the hadrons. Recently, a complete TMD- 
factorization derivation, in terms of well-defined TMDs 
with individual evolution properties, was presented by 
Collins in Ref. 14[. Refs. Ia,ll6| further illustrated how 
the formalism can be used to obtain evolved TMDs from 
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fixed-scale fits for unpolarized TMDs and the Sivers func- 
tion. 

Existing extractions of the Sivers function using the 



TSSA, A 



sin(0 h — 0s) 
UT 



have been performed using experi- 
17H22J. These extractions 



mental data at fixed scales 
provide interesting information about TMD effects at the 
fixed scales where they are performed; however, without 
a reliable way to evolve them to different scales, their 
predictive power remains limited. 

The purpose of this letter is to demonstrate that by 
using QCD evolved TMDs one can explain an observed 
discrepancy between HERMES and COMPASS data and 
for the first time make predictions for upcoming experi- 
ments at higher energy scales on the basis of a complete 
and correct treatment of evolution for the TMDs. 

Definitions and Notation: The differential cross 
section for SIDIS, 1(1) + N(P,S) -> 1(1') + h{P h ) + X 



is 



da 



dxdydzdtfihdifisPh-LdP) 
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ML^W^ (1) 



where Phj_ is the transverse momentum of the final state 
hadron h, and where we utilize the standard kinematical 
variables: q 2 = -Q 2 , x = Q 2 /2P ■ q, y = P ■ q/P ■ l,z = 
P ■ P h /P ■ q. The TMD-factorization formula for SIDIS 
in terms of well-defined TMD PDFs is E 



f 

x J d 2 p T d 2 K T 5 {2) (zp T + K T - P h ±) 

x F f/Pt (x, zp T ,S;p, ( F )D h/f (z, K T ; p, ( D ) 

+ Y(P h± ,Q), (2) 

where all non-perturbative information is encoded in 



the TMD PDF F 



f/Pt 



and the TMD FF D h/f while 



\Hf(Q 2 ,fJ,)\ ,lu is a perturbatively calculable hard part. 
The Y(Ph±, Q)-term gives the correct treatment of the 
cross section at high Ph± ~ Q in terms of collinear factor- 
ization. As is common, the renormalization scale is set to 
/i = Q. The parameters £f,Cd, which are related to the 
rcgularization of rapidity divergences, obey CfCd ~ Q 4 - 
(Consistency with perturbation theory requires p, a/Cf 
and \/Cd each to be of order 0(Q).) 

The Sivers asymmetry is defined as a the ratio of cross 
section combinations: 

,sin(4> h -4>s) _ 

J d<j> h d(j) s 2 sm(<j) h - <l>s)(cr(<t>h,4>s) ~ &((f>h,4>s + tQ) 



/ d<t> h d<t>s(<j((t>h,<l>s) + <*(4>h,<t>s + i")) 



(3) 



In the numerator, the integration over azimuthal angles 
with a sin(<j>h — 4>s) weighting factor projects out the 



Sivers effect. The numerator and denominator may also 
be integrated over x, z and/or Ph± depending on the 
particular combination of variables one is interested in. 

The asymmetry is obtained by applying the TMD- 
factorization in Eq. ^ to obtain cross sections in Eq. ([3]) . 
The calculations themselves are typically done in trans- 
verse coordinate &T-space in terms of structure functions, 
whose relations to the differential cross section are given 
in Ref. [26[ • In the case of the Sivers function, the general 
expression for the evolved TMD in coordinate space was 
found in Ref. [l6| to be 

F£ f (x,b T ;Q,( F ) = F^ f (x,b T ;Q 0l Q 2 ) 

~ dp,' 
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Analogous formulas hold for the unpolarized TMDs. 
The symbols, jk and j F , are perturbatively calculable 
anomalous dimensions, K(b*;fib) is the perturbatively 
calculable CS kernel written in terms of 6* which is the 
prescription for matching to the region where 1 /bp can be 
treated as a perturbatively large scale. We use the usual 
prescription of 0] where 6* = bp/ \J\ 



b 2 T /bl 



and 



fib = C\/b*, and b max and C\ are parameters to be spec- 
ified later. Note that it is the derivative F^, (x, bp; Q, £) 
of the QCD evolved coordinate-space Sivers function with 
respect to bp that appears in Eq. ([4]) for the evolution. 
Qo is the starting scale for the evolution. The non- 
perturbative but universal and scale-independent func- 
tion gxipp) describes the behavior of K(bp\pb) in the 
non pcrturbative region at large bp. An important 
prediction from the TMD-factorization theorem is that 
gxibp) is universal, not only between different processes, 
but also between all different types of quark TMDs (both 
PDFs and FFs). 

For this letter, we assume that Q is low enough that 
we can neglect the F-term in Eq. ([2]) [l5|. Furthermore, 
we use a Gaussian ansatz to parametrize the input distri- 
bution F^t ^(x, bp; Qo, Qq), though this means that we 
do not utilize the fact that at larger P/ l ± ^> Aqcd the 
TMD PDFs are related to collinear distributions through 
perturbative coefficient functions. (In the Sivers case, 
this involves the Qiu-Sterman function [HEl.) Still, m 
Ref. [15^ it was shown that a Gaussian ansatz provides 
a good description of the evolved Sivers function for the 
low region of transverse momentum and moderate hard 
scales we are interested in for this letter. Several groups 
have parametrized the polarized and unpolarized TMD 
PDFs and FFs at fixed scales in terms of simple Gaussian 
HESS!, 

l31j . and these may be used as the input 



fits 



functions for the evolution. 
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FIG. 1: Comparison between HERMES '27] and preliminary COMPASS data fH for the (a.) z and (b.) P h± dependence of 
Eq. © with a proton target and tt + and h + as final state hadrons respectively. The solid line is the fit from Ref. [53] . The 
dashed curve is the result of evolving to the COMPASS scale using the full TMD-evolution of Ref. 



Analysis and Discussion: As input distributions, 
we use the already existing Gaussian paramctrizations of 
the Torino group 0, relevant for low (Q 2 )Hcrmcs — 2.4 
GeV 2 and typical for the HERMES experiment. These 
earlier fixed scale fits were done at leading order in QCD 
and neglecting the QCD evolution of the TMDs, which 
was not available at that time. We note that the anal- 
ysis of Ref. [22| also uses deuteron data §52j from the 
COMPASS experiment, which corresponds to higher val- 
ues of Q 2 . However, the COMPASS asymmetry [32| on 
the deuteron target is very small due to strong cancel- 
lations between the up and down quark Sivers functions 
and thus is not heavily affected by the evolution. We 
have verified that the results of the Torino fits are negli- 
gibly altered if the deuterium data are excluded and only 
HERMES data [ItJ are used in the fit, and the main re- 
sult of our present analysis is not affected. 

Our calculations will follow the steps of Ref. . For 
gK, we use the functional form gx = \g2^\ with 172 = 



0.68 GeV 2 [33j , which was obtained by fits performed 
using Drell-Yan data. In Eq. ([4]), this corresponds to 
using C\ = 1.123 and b max = 0.5 GeV -1 . All anomalous 
dimensions and K are calculated to lowest non- vanishing 
order as in Refs. (l4, 15|. 



In Fig. nja,b), we show the evolution using the full 
TMD-factorization approach as expressed in Eq. ([4]), 
where the evolution is due to the terms in the expo- 
nential. The evolution is applied to the most recent 
Torino fits 22j as a function z and Ph± , and use 



hard scales corresponding to both HERMES data [27 1 
and recent preliminary COMPASS data jH]. The re- 
sult of the evolution is compared with the data. The 
x-dependent asymmetry is not ideal for the comparison 
because there are strong correlations between x and Q 2 . 
(Recall Q 2 ~ xys.) However, z or P^± dependent asym- 
metries are measured at almost the same hard scales, 



namely (Q 2 )Hormos ^ 2.4 GeV 2 and (Q 2 )compass - 3.8 
GeV 2 , so we focus on the Sivers asymmetry as a func- 
tion of these variables. (For the preliminary h + COM- 
PASS data that we use, (Q 2 ) varies between 3.63 GeV 2 
and 3.88 GeV 2 , in the range of z from 0.2 to 0.7. The 
corresponding variation in our calculation is negligible 
relative to the variation between the HERMES and pre- 
liminary COMPASS data sets.) We observe that includ- 
ing QCD evolution leads to excellent consistency between 
the HERMES ^ and preliminary COMPASS data [H], 
without the need for further fitting. The two data sets 
correspond to different ranges in x, and this could be 
partly responsible for the variation. A similarly fast evo- 
lution has not been seen so far in the Collins Single Spin 
Asymmetry 2^, 34|, suggesting a more complicated in- 
terplay between bx, x and z dependence. We leave a 
careful consideration of these issues to future studies. 
Nevertheless, we find the early success of the compari- 
son in Fig . (fTI) encouraging, especially as leading order 
fits [13, wL [H] fail to reproduce COMPASS proton data 
28] sufficiently well. Still, we caution that future fits will 
need to account for the x-dependence as well. 

A critical point is that the information about the non- 
perturbative evolution contained in gx is taken from the 
measurement of a totally different observable, at much 
higher energy scales [33[ (unpolarized Drell-Yan scatter- 
ing up to Tevatron energies). In Fig. QJb) we show a 
similar plot but for the Ph± dependence. That the same 
g K successfully describes TSSA at HERMES and COM- 
PASS is compelling evidence for the universality of gx 
predicted by the TMD-factorization theorem. 

In Fig. [2l we show the evolution of the full asymmetry 
to higher values of Q 2 . Note that, although Refs. [15l.ll6| 
report a strong suppression of the unpolarized TMDs and 
the Sivers function itself with increasing Q 2 , the TSSA is 
not as heavily suppressed. Therefore, it may be expected 
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FIG. 2: Sivers evolution in Q , integrated over x, z and Pha 



that the Sivers SSA remains significant at the higher Q 
values of experiments planned at the Rclativistic Heavy 
Ion Collider (RHIC) and the EIC. Still, the QCD evolu- 
tion effects are clearly non-negligible and should be cor- 
rectly included in future extractions. Ref. Q predicts a 
roughly ~ 1/yfQ suppression for the peak of the Sivers 
asymmetry as a function of transverse momentum, for 
large Q 2 > 10 GeV 2 . We find that, for the full asymme- 
try integrated over all transverse momentum, a power- 
like scaling law does not provide a good description in 
the range of Q in Fig. [U Generally, we find that the evo- 
lution leads to suppression that is faster than ~ 1/y/Q, 
but slower than ~ 1/Q 2 - We caution, however, that a 
completely correct treatment at large Q must include the 
Y-term in Eq. , and it is possible that this will weaken 
the rate of the suppression. 

To conclude, we remark that it is important for future 
theoretical calculations to not only explain experimen- 
tal results, but also to make precise pQCD-based pre- 
dictions that can be tested against future data at larger 
Q. With this in mind, we view the success of the TMD- 
factorization treatment in explaining the HERMES and 
COMPASS as highly encouraging. 
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